Time Evolution of the Fine Structure Constant in a Two-Field Quintessence Model 
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We examine the variation of the fine structure constant in the context of a two-field quintessence 
model. We find that, for solutions that lead to a transient late period of accelerated expansion, it 
is possible to fit the data arising from quasar spectra and comply with the bounds on the variation 
of a from the Oklo reactor, meteorite analysis, atomic clock measurements. Cosmic Microwave 
Background Radiation and Big Bang Nucleosynthesis. That is more difficult if we consider solutions 
corresponding to a late period of permanent accelerated expansion. 
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I. INTRODUCTION 



The recent claim that the spectra of quasars (QSOs) 
indicates the variation of the fine structure constant, a, 
on cosmologically recent times [J 0, S has raised con- 
siderable interest in examining putative sources of this 
variation. In most models a possible variation of the 
fine structure constant is studied by arbitrarily coupling 
fields to electromagnetism, as suggested by Bekenstcin 
0- Thus the proposals put forward sofar consider a 
scalar field Q (with an additional coupling to dark 
matter and to a cosmological constant in the former case) 
and quintessence 's', In fact, as discussed in Ref. [loj . 
the couplings of gravito-scalar fields, such as the axion or 
the dilaton, to electromagnetism naturally arise in = 4 
Supergravity in four dimensions, making this model par- 
ticularly interesting. It is worth mentioning that, in the 
latter model, the mass of these scalars can drive the ac- 
celerated expansion of the Universe, transient or eternal 

m 

In this work, we shall consider the implications for the 
variation of a of a two-field quintessence model (l j| , with 
the quintessence fields coupled to electromagnetism, as 
proposed by Bekenstein There are several motiva- 
tions for studying potentials with coupled scalar fields. 
Firstly, if one envisages to describe the Universe dynam- 
ics from fundamental theories, it is most likely that an en- 
semble of scalar fields (moduli, axions, chiral superfields, 
etc) will emerge, for instance, from the compactification 
process in string or braneworld scenarios. Furthermore, 
coupled scalar fields are invoked for various desirable fea- 
tures they exhibit, as in the so-called hybrid inflationary 
and reheating models PI. The model of. Ref. 12] has 
the additional bonus of leading to transient as well as 
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permanent solutions for the late time acceleration of the 
Universe. The former solutions are desirable given that 
it has been recently pointed out that an eternally acceler- 
ating universe poses a challenge for string theory, at least 
in its present formulation, since asymptotic states are in- 
consistent with spacetimes that exhibit event horizons 
p^ . Moreover, it is argued that theories with a stable 
supersymmetric vacuum cannot relax into a zero-energy 
ground state if the accelerating dynamics is guided by 
a single scalar field .141] . a problem that can be circum- 
vented in the two-field model we are considering p^ . 

We now turn to the available observational bounds on 
the variation of a. Observations of the spectra of 128 
QSOs with z = 0.2 — 3.7 suggest that, for z > 1, a was 
smaller than at present P, [ 



Aa _ a{z) 
a 



ao 



= (-0.54±0.12) X 10 



-5 



(1) 



at A.7a. 

The most recent data is from Chand et al. 0, 0| 
obtained via a new sample of Mg II systems from distant 
quasars with redshifts in the range 0.4 < z < 2.3 yield 
more stringent bounds (3 cr), namely: 



Aa 

a 



(-0.06 ±0.06) X 10' 



(2) 



where terrestrial isotopic abundances have been assumed. 
If, instead, low-metalicity isotopic abundances are as- 
sumed, Chand et al. obtain 



Aa 



(-0.36 ±0.06) X 10" 



(3) 
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in which case the statistical inconsistency with Murphy 
et al., Eq. is clearly smaller. Notice that, in contrast 
with Webb et al., who use different lines from different 
mtiltiplets and elements, Chand et al. use mostly Mg II 
data yielding a smaller but better quality dataset. 

On the other hand, the Oklo natural reactor provides 
a bound, at 95% CL, 



2 



-0.9xl0"^< < 1.2 X 10"^ , (4) 

a 

for ^ = 0.14 'SHI m. Notice, however, that the use of 
an equilibrium neutron spectrum has been criticized in 
Ref. |20| , where a lower bound on the variation of a over 
the last two billion years is given by 



ai z = 10^ [23, , and from Big Bang nucleosynthesis 
(BBN), 

- 6 X 10""^ < Aa/a < 1.5 X 10"^ , (11) 

at z = IQS - 10^° IISIS^- 

In addition, we should take into account the equiva- 
lence principle experiments, which imply Q 



Aa 



> +4.5 X 10~ 



(5) 



Estimates of the age of iron meteorites (2 = 0.45), 
combined with a measurement of the Os/Re ratio result- 
ing from the radioactive decay ^^^Re ^*^0s, gives 



at la, and 



(-8 ±8) X IQ-^ 



7 Aa 

- 24 X 10~^ < < 8 X 10 

a 



-1-7 



(6) 



(7) 



at 2(7. 

Notice that, if the variation of the fine structure con- 
stant is linear, the Murphy et al. observations of QSO 
absorption spectra and of geochemical tests from me- 
teorites are incompatible given that the former yields 
Aa/a ~ 5 X 10~^ at z = 0.5, while the latter leads to 
Aa/a < 3 X 10"^ ai z — 0.45 However, this problem 
may not exist if one uses the Chand et al. dataset. 

Moreover, observations of the hyperfine frequencies of 
the ^"^^Cs and ^^Rb atoms in their electronic ground 
state, using several laser cooled atomic fountain clocks, 
give, at present (z — 0) (see also psf') 



< 4.2 X 10"^^ 



yr 



(8) 



where the dot represents differentiation with respect to 
cosmic time. Tigher bounds arise from the remeasure- 
ment of the Is — 2s transition of the atomic hydrogen and 
comparison with a previous measurement with respect to 
the ground state hyperfine splitting in ^■^^Cs and combi- 
nation with the drift of an optical transition frequency in 
i99Hg+, that is 



- = (-0.9 ±4.2) X lO"^'^ yr~i 
a 



(9) 



There are also constraints coming from Cosmic Mi- 
crowave Background Radiation (CMBR) , where 



|Aa/a| < 10" 



(10) 



< 10" 



(12) 



where C,p is the coupling between the scalar and elec- 
tromagnetic fields and u = ^ip"; ^ being the reduced 
Planck mass (M = Mp j VStt) and A/» the corresponding 
analogue in the scalar sector. 

Notice that not all models put forward sofar satisfy 
the abovementioned bounds. For instance, quintessence 
models like the last one in ^ and the = 4 Supergravity 
models in four dimensions are not consistent with Webb 
et al. data. In the following, we shall present our model 
and show that it is consistent with all the available data 
for a suitable choice of the model parameters. 



II. MINIMAL COUPLING OF QUINTESSENCE 
TO THE ELECTROMAGNETIC FIELD 

We shall study the time evolution of a in a model where 
two homogeneous scalar fields are minimally coupled to 
electromagnetism. This evolution follows from the effec- 
tive action, in natural units {M = 1), 



S 



--R + Lb + Cq + C. 



(13) 



where Cb is the Lagrangian density for the background 
matter (CDM, baryons and radiation), which we con- 
sider homogeneous, corresponding to a perfect fluid with 
barotropic equation of state pb = Wb Pb, with Wb constant 
(—1 < Wb < 1; Wb = 1/3 for radiation and Wb = for 
dust); Cq is the Lagrangian density for the scalar fields 



We consider the two-field quintessence model with po- 
tential m 



y(0,V) = e-^M0» 



(15) 



where 



P((/), tP) = A + {(!)- 0,)2 + B ii;- 

+ C (t){ij- + Di}{(j)- (j>,f . (16) 
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FIG. 1: Evolution of the quintessence fields (upper panel), 
the equation of state parameter (middle panel) and the 
quintessence fractional energy density (lower panel) , for tran- 
sient (dashed) and permanent (full) acceleration solutions. 

Potentials of the type exponential times a polynomial, 
involving one single scalar field, have been considered be- 
fore, see e.g. lU- The evolution equations for a spatially- 
flat Friedmann-Robertson- Walker (FRW) Universe, with 
Hubble parameter H = a/a, are 

H = (pb+Pfc + <^'+^') , 

^ = -ZH^-d^V , (17) 
subject to the Friedmann constraint 

= \{p^ + \4>' + + , (18) 

where V = • The total energy density of the 

homogeneous scalar fields is given by pg = 0^/2 -I- -0^/2-1- 
!/((/., V). 

Integrating Eqs. IjlTI) . one finds that there are essen- 
tially two realistic types of solutions, corresponding to 
either transient or permanent acceleration, as illustrated 
in Figure n It is possible to obtain permanent or tran- 
sient vacuum domination, satisfying present bounds on 
observable cosmological parameters, for a rather broad 
range of parameters of the potential, as discussed in ■ 
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FIG. 2: Evolution of a for a transient acceleration model 
with = 2 X 10"*^ and C2 = 8 x 10~^. In the upper panel, 
the boxes represent the QSO bounds given in Eqs. ^ (top 
box), (middle box) and Q (lower box). Also shown is 
the QSO absorption systems dataset of Chand et al. (Table 
3 of Ref. tl5|)- The lower panel details the behaviour of 
a for small values of z. The lower left plot shows the Oklo 
bound (dash-dotted lines) and the right one the meteorite 
bounds (dash-dotted and dotted lines correspond to la and 
2(T, respectively). 

Permanent vacuum domination takes place when at 
least the 4> field settles at the minimum of the poten- 
tial; soon after, ?/; monotonically increases until it reaches 
its asymptotic value. For transient vacuum domination 
(which occurs either when the potential has no local min- 
imum or when arrives at the minimum with enough 
kinetic energy to roll over the potential barrier) the 
field increases monotonically and -0 slightly decreases. 
The contribution of P(0, tp), which controls the presence 
and location of the minimum, is only important at re- 
cent times, since at early epochs the exponential factor 
dominates and the model behavior depends essentially 
on parameter A. Notice that, as A increases, the model 
becomes more similar to ACDM at early times. 

The interaction term between the scalar fields and the 
electromagnetic field is, as suggested in Ref. 0, 

-Cem - -\Bf{4>, V^)Fm-F'" > (19) 

with a linearly expanded i?(0, i/") given by 

(0, 0) = 1 - Ci(0 - 0o) ~U^- ^0), (20) 

where 00 ■00 ^-re the present values of the scalar fields. 
Therefore, the variation of the fine structure constant, 
a = ao/Bpi'p^i'), is given by 



4 




0.05 0.1 0.1 0.2 0.3 0.4 

z z 



FIG. 3: As for Fig.|5|but for a permanent acceleration model 
with = -4 X 10"^ and C2 = 1 x 10"^ 



— =Ci(0-0o) + C2(V'-^o). (21) 
a 

Notice that the upper bound on Cf, see Eq. (I12II . im- 
phes that Ci,2 < 7 x 10~^ because of the different defini- 
tion for Bp used in 0, namely, Bf{^) = 1 — -^{ip— ipo). 

The rate of variation of a, at present, can be written 

as 

where y = 1 + z and Hq is the value of the Hubble con- 
stant at present, Hq = {h/9.78) x 10^^ yr"^ 



III. RESULTS 

We adopt the method of choosing the model param- 
eters so as to satisfy a set of priors for the cosmologi- 
cal parameters {h — 0.70, — 0.3, Slg — 0.70 and 
= 4.15 X lO^^h^^) and then adjust the coupling pa- 
rameters, Ci and C2, in order to satisfy the bounds on 
the evolution of a. Notice that the priors chosen above 
are consistent with a combination of WMAP data and 
other CMB experiments (ACBAR and CBI), 2dFGRS 
measurements and Lyman a forest data, which gives 32] : 
h = 0.7ltj]:[Jt = 0.27 ± 0.04, riQ = 0.73 ± 0.04 and 
WQ < —0.78 (95% CL), corresponding to the best fit for 
the observed Universe. 
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FIG. 4; Evolution of a for a transient acceleration model with 
Ci = 2 X 10"'' and C2 = 8 X 10"^ (fuU line) , d = 5.3 x 10"^ and 
C2 = 3 X 10"^ (dashed line), Ci = 1-4 x 10~^ and C2 = 7 x 10~* 
(dash-dotted line). Line and box conventions are those of 
Figs. Hand 01 



For large z, the tightest bound on dark energy arises 
from nucleosynthesis, Jig (2 = lO^*') < 0.045, implying 
that A > 9 l3fli |. Considering a possible underestimation 
of systematic errors, one gets a more conservative bound, 
Hq < 0.09 at 2(7, which corresponds to A > 6.5. For the 
model parameters chosen in this work (see below) , we get 
fig = 0.042. 

The model is also consistent with CMBR observations. 
Recent CMBR data imply that Hq < 0.39 at last scatter- 
ing _30jj which is less stringent than the nucleosynthesis 
bound and is clearly obeyed by our model. We have 
also computed the angle subtended by the sound hori- 
zon at last scattering, obtaining Oa — 0.599°, clearly 
within WMAP bounds on this quantity 32]: 6a = 
0.598 ± 0.002°. 

As an example of transient acceleration models, we 
consider the following set of parameters: A — 9.5, A = 
0.1, B = 10-3, C = 8 X 10-^ D = 2.8, 0* = 28.965, 
-0* = 20, with Ci = 2 X 10~s and C2 = 8 x 10"^ In this 
case, the variation of a is in agreement with both the 
Oklo and meteorite bounds, as can be seen in the lower 
panel of Figure |2| In the upper panel, we see that the 
fit to Chand et al. QSO data (where terrestrial isotopic 
abundances have been assumed) is acceptable: we get 

/d.o.f = 22/15 (notice that, in the computation of the 
d.o.f., only 8 free variables should be taken into account 
since one potential parameter is tuned in order to have 
the observed dark energy density at the present). Also 
the constraint from atomic clock measurements is verified 
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It is clear that the evolution of a, past and future, is 
a reflection of the dynamics of the fields (j) and ip, which 
depends crucially on whether acceleration is permanent 
or transient. In the example of Figure |21 a is increas- 
ing monotonically and will continue to do so. However if 
Ci ^ C21 it is the ip field that determines the evolution; 
hence, in that case, a may decrease in the future. For per- 
manent acceleration solutions, the sign of (^2 determines 
the future evolution of a, since ip has more dynamics 
than (p. Hence, for the example of Figure |31 a will in- 
crease in the future since (2 > 0; the oscillating behavior 
of a reflects the oscillations of cj) before it settles at the 
minimum of the potential. 

In Figs. ElandElwe show that, increasing the couplings 
to the electromagnetic field, and (2, it is also possible 
to fit Murphy et al. QSO data [H, both in the transient 
and permanent acceleration regimes. However, as should 
be expected, it is more difficult to respect the Oklo and 
meteorites bounds in both cases. 



FIG. 5; As for Fig.0]but for a permanent acceleration model 
with (^1 = -4 X 10"^ and C2 = 1 x 10"® (full line), Ci = 
-1 X 10"* and C2 = 1 x 10"® (dashed line), Ci = -1-5 x 10"* 
and C2 = 1-4 x 10"^ (dash-dotted line). 

as we obtain 



- = -4.5 x 10-1^ yr-i . (23) 
a 

Moreover, the model gives variations in a during BBN 
and at the last scattering surface well within the bounds 
stated earlier, e.g. we get 



(CMBR) = -2.7 x 10"^ , 

a 

— (BBN) = -1.1 X 10"^ . (24) 
a 

For models with permanent acceleration, however, con- 
sistency with QSO spectra and Oklo data is problematic. 
Our results arc illustrated in Figure 13 where we have 
chosen A = 9.5, A = 0.02, B = 2x IQ-^, C = 6 x IQ-^^ 
D = 4.5, 0, = 28.9675, -0, = 15, for Ci = -4 x 10"''^ and 
C2 = 1 X 10"^. We get 

— = 5.2 X 10"^^ yr~^ 
a 

— (CMBR) = 4.5 X 10"'^ 
a 

— (BBN) =2.9x10-^ (25) 
a 

and x^/d-o.f = 21/15. Notice that the BBN bound is 
not respected in this case. 



IV. CONCLUSIONS 

In this work, we have studied the variation of the fine 
structure constant in the context of a quintessence model 
with two coupled scalar fields. We find that transient ac- 
celeration models can fit the latest QSO data and com- 
ply with the upper bounds on Aa from the Oklo reactor, 
meteorite analysis and the atomic clock measurements. 
For permanent acceleration models, however, it is more 
difficult to fit the QSO data and satisfy the Oklo, me- 
teorite and BBN bounds simultaneously. We have stud- 
ied the sensitivity of our results to Ci and (2, the cou- 
plings of quintessence fields with electromagnetism, and 
we have found that, in order to be consistent with the 
data, these parameters must be at least one order of 
magnitude smaller than the upper bound implied by the 
Equivalence Principle. 

On a more general ground, we could say that estab- 
lishing whether there is a variation of the fine structure 
constant and, in the affirmative case, identifying its ori- 
gin, remains a difficult task before a deeper analysis of 
the systematic errors of the observations and studies of 
the degeneracies with the various cosmological param- 
eters. This is particularly evident in what concerns the 
compatibility of the datasets of Murphy et al. and Chand 
et al. 
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